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CN ; ABSTRACT 

^T^ Aims. With the first ~ 10000 spectra of the flux limited zCOSMOS sample {Iab ^ 22.5) we want to study the evolution of environmental effects on 

^^ galaxy properties since z ~ 1.0, and to disentangle the dependence among galaxy colour, stellar mass and local density. 

r \ Methods. We use our previously derived 3D local density contrast 6, computed with the 5'* nearest neighbour approach, to study the evolution 
• , with z of the environmental effects on galaxy U-B colour, D4000A break and [011]/i3727 equivalent width (£^[011]). We also analyze the impu- 
re cations due to the use of different galaxy selections, using luminosity or stellar mass, and we disentangle the relations among colour, stellar mass 
C^ and 6 studying the colour-density relation in narrow mass bins. 

JL Results. We confirm that within a luminosity-limited sample (Mg < -20.5 - z) the fraction of red (U - B > 1) galaxies depends on S at least 

i , up to z ~ 1, with red galaxies residing mainly in high densities. This trend becomes weaker for increasing redshifts, and it is mirrored by the 

"*T* behaviour of the fraction of galaxies with D4000A break > 1.4. We also find that up to z ~ 1 the fraction of galaxies with log(£lV[0//]) > 1.15 is 

^ higher for lower 6, and also this dependence weakens for increasing z- Given the triple dependence among galaxy colours, stellar mass and S, the 

I I colour-(5 relation that we find in the luminosity-selected sample can be due to the broad range of stellar masses embedded in the sample. Thus, we 

study the colour-5 relation in narrow mass bins within mass complete subsamples, defining red galaxies with a colour threshold roughly parallel 

T-H to the red sequence in the colour-mass plane. We find that once mass is fixed the colour-5 relation is globally flat up to z ~ 1 for galaxies with 

^ log(M/MQ) > 10.7. This means that for these masses any colour-(5 relation found within a luminosity-selected sample is the result of the combined 

T-H colour-mass and mass-^ relations. On the contrary, even at fixed mass we observe that within 0. 1 < z < 0.5 the fraction of red galaxies with 

^^ log(M/Mo) < 10.7 depends on S. For these mass and redshift ranges, environment affects directly also galaxy colours. 

OO Conclusions. We suggest a scenario in which the colour depends primarily on stellar mass, but for an intermediate mass regime (10.2 < 

C^ log(M/Mo) < 10.7) the local density modulates this dependence. These relatively low mass galaxies formed more recently, in an epoch when 

I * more evolved structures were already in place, and their longer SFH allowed environment-driven physical processes to operate during longer 

L«. periods of time. 

^^ Key words. Galaxies: evolution - Galaxies: fundamental parameters - Galaxies: statistics - Galaxies: high-redshift - Cosmology: observations - 
T-H . Large-scale structure of Universe 
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1. Introduction properties have been used in order to classify galaxies. The 

inter-correlations and scaling laws that were found among these 
Galaxies can be considered the building blocks of the universe, properties and also with morpholo gy (to give only a fe w ex- 
and nevertheless we are far from a complete understanding of ^j^ples, the 'Fundamental Plane', iDiorgovski & D^ [19^ 



their formation and evolution. Although theoretical scenarios Dressier ^tiD SMI the 'Photometric Cube', ' Scodeggio et al. 
can help us to address these open issues, only through observa- oqq^. the colour- magnitud e relation. Wi^^^alhm & Sandal 
tions we can robustly shape a coherent picture, following back- [iqttI [TuUv et al.l i 1982; GavazzietalJ [19961) have heelTmn- 



ward galaxy lifetimes from what we seem the local universe up ^^^^^ j„ j^e last years using the data provided by large 

to their undisclosed origins. Jq^.^1 galaxy surveys like the Sloan Digital Sky Survey 

The first observations of near by galaxies led to t heir mor- (SDSS iYork etaP llOOOh and the Two degr ee Field Galaxy 

phological classification (|Hubble| [1926, but see also [SandaM Redshift Survey (2dFGRS,'Colless et al."2001'), with the studies 

[19751). In more recent years, several internal physical galaxy bijKauffmann et al. (2004), Baldry et al. (2004), Balog h^TZI 

(12004 ) and several others. These findings suggested a broad bi- 

Send offprint requests to: O. Cucciati (olga.cucciati@oamp.fr) modal galaxy classification: 1) the red, passive galaxies, with el- 

* European Southern Observatory (ESO), Large Program 175.A- liptical/spheroidal shapes and composed by an old (on the mean) 
0839 
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population of stars, moving mostly on disordered radial orbits; 
2) the blue, active galaxies, disk-dominated, where stars (of rel- 
atively young age) are concentrated in circular orbits. 

The close relationship between these galaxy types and lo- 
cal environment was primarily uncovered via the study of 
nearby clusters. T he so-cal l ed m orpholo gy-dens i ty rela tion, as 
first quantified by lOemled (119741) and IDressled (ll980), holds 
that star-forming, disk-dominated galaxies tend to reside in re- 
gions of lower galaxy density, such as cluster outskirts and the 
field, while red, elliptical galaxies are found preferentially in 
the higher density regions, like cluster cores. Recent work based 
on the 2dFGRS and the SDSS has established that the connec- 
tions between local environment and morphology holds not only 
in clusters, but extend continuously over the full range of local 
densities, from the centers of clusters out into the field popula- 
tion. Moreover, not only morphology, but several other galaxy 
properties correlate with environment. On average, high den- 
sity regions are populated mostly by redder, brighter and less 
star forming galaxies, w hile the opposite is true for the low 
density environment (e.g. Balogh et al. '2004"; 'Kauffmann et alJ 
[2004; Tanaka et al. 2004; Blanton et al. 2005; Hogg et al. 2004). 
These studies also indicated that colour is the property most 
tightly related to environment. 

With the recent advent of large deep galaxy surveys, it has 
been found that the galaxy bimodal distribution is not charac- 
teristic of the local universe only. On the one side, early type 
galaxies are found to be the predominant population among red 
and lumi nous galaxies at all red shift up to z ~ 1.3 (see e.g. 
iBell et a l. 2004; Menan teau et al . 2006). On the other side, it has 
been observed that the bimodal colour distribution holds even 
up to z ~ 2, although with a significant contrib ution of dusty 
star forming galaxies to th e red population (e.g. . IWillmer et alJ 
120061; iFranzetti et al.ll2007l) . 

Also the environment has been shown to extend its influence 
on galaxy properties at these high redshifts. In the last decade 
studies of increasingly distant clusters (to z ~ 1) have grown 
in number, thanks to high-resolution imaging and spectroscopic 
data. These works have shown, for example, that the fraction of 
early type galaxies d ecreases fr om the inner part of clusters to 
their outskirts (e.g. . iTreu et an i2003). and the fraction of star 
forming ga laxies is higher in the field than in any other environ- 
ment (e.g. . iPoggianti et al.ll200d) . Moreover, it has been shown 
that, at least up to z ~ 1, the fraction of blue galaxi es anti- 
corre lates with group richness (lo vino et al. 2010 ; Cucciati et alJ 
1201 0|). while the fracti on of early-type galaxies correlates with 
it ('Kov ac et al.ll2010bl) . Still, comparisons of local results with 
studies of high-redshift clusters have pointed towards significant 
evolution in the relationship between galaxy properties and lo- 
cal environment from z ~ 1 to the present. For example, the 
morphology-density relation for galaxies i n clusters and group s 
evolves strongly with increasing redshift dPressler et al.|[l997l) . 
and there is evidence that the evolution in intermediate density 
environments has occurred more recently than in high density 
environments (ISmith et al.l 12005.) . It also appears that the mor- 
phological mix of cluster early type galaxies (the relative frac- 
tion of le nticular and elliptical galaxies) has changed with time 
dFasano e t al. 2000; Postman et al. 2005). 

Nevertheless, with these high redshift cluster studies it 
is possible to divide galaxies only into broad environment 
classifications (such as field, group and cluster populations). 
Furthermore, clusters include only a relatively small fraction (by 
number and an even smaller one by volume) of the total galaxy 
population at any epoch, and therefore they offer only a limited 
vision on galaxy evolution. The new generation of large and deep 



redshift surveys very recently allowed us to extend to high red- 
shift the investigation of environmental effects on several galaxy 
properties within the full range of local densities, from voids 
and very low densities up to the higher density peaks in the cen- 
ter of clusters. These high redshift studies could help us resolve 
the still open debate about the origin of environmental effects 
on galaxy properties. The question is whether these properties 
were imprinted upon the galaxy population at their formation 
epoch (the so called 'nature' hypothesis), or whether they are a 
result of environment-driven evolution, i.e. the end products of 
processes that have operated over a Hubble time (the so called 
'nurture' scenario). 

Many physical mechanisms that could be responsible for 
this correlation between galaxy properties and environment 
have been proposed, but up to now there is no a clear un- 
derstanding of their real role on the observed trends. From 
preliminary studies of local density and group membership at 
Z ~ 1 that have been developed independently and nearly at 
the same time in the se past few years b y the VIMOS-VLT 
Deep Survey (VVDS, Le Fevre et al. 2005) and by the DEEP2 
Galaxy Redshift Survey (.Davis et al.. .2003.) . it seems at least 
that both a biased galaxy formation (galaxies formed earlier in 
high matter-density peaks, giving galaxies different imprinting 
of initial conditions, Marinoni et al. 2005) and a complex evo- 
lution of star formation quenching (for example the weakening 
of the colour-dens ity relation at higher redshift, ICucciati et al.l 
|2006 | and Cooper et al. 20071 with possibl y its reversal at z 



1.5, Cucciati et al. 2006, but see also Elbaz et aP l2007l and 
Cooper et al. 2008) concurred to build up the observed relations 
between environment and galaxy properties. We will discuss 
more deeply high redshift results in the following of this paper 

A detailed view of the environmental effects on galaxy prop- 
erties is still missing, however. The scale length at which en- 
vironment affects the most galaxy properties has been studied 
in the local universe, and it seems that local density acts on 
colour. H a equivalent width and D4000A b reak on a s cale of ~ 1 
Mpc /i ' (iKauffmann et al.l2004 but also .Blanton et al..2006. for 
scales also ~ 0.5 Mpc h'^), and what it is observed on larger 
scales is only a mirror of what really happens on smaller dis- 
tances. Nevertheless, these results are related only intuitively at 
the typical scales of cluster halos, and a precise quantification 
of these observed effects has still to be done. Moreover, as it is 
known that several properties are inter-related and they also are 
related with environment, one would like to understand whether 
the local density affects only a property, that then drives the 
dependence of the others, or more than one. For example, we 
already know that mass is responsible for shaping gala xy star 
formation history jG avazzi & Scodeggio 1996; Scodeggi o et al.l 
I2OO2I; IKauffmann et al.ll2003l) . and moreover mass does depend 
on environment, as it as been shown at b oth local and high 
(z ~ 1) redshi ft (Kauffmann et al. 2004 ijBundv et al.1 120061: 
IScodeggio etani2b09t iBolzonella et al.ii2009h . Should we con- 
sider the galaxy stellar mass as the only property directly af- 
fected by the environment, this meaning that the other prop- 
erties depend on local density only through their dependence 
on mass? The answer seems not to be unique, as at z ~ 0.1 
the galaxy colour, D4000A break and sSFR are found to de- 
pend oneimronm^nteven when mass is fixed (Kauffmann et al] 
l2004t iBaldrv et al.l l2006h . but no colour-density relation has 
been found in the range 0.2 < z < 1 .4 in mass bins within VVDS 
data dScodeggio et al.l2009l) . 

In this frame of open issues, we use the first ~ 10000 
zCOSMOS spectra to inspect more deeply the environmental 
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effects up to z ~ 1. With its pure flux limit {Iab < 22.5) and 
small redshift measurement error (~ lOOkm/s), this galaxy sam- 
ple allowed us to reliably reconstruct the local density field o n 
small scales up to z ~ 1, as described in iKovac et alJ (1201 Oah . 
Moreover the quality of the spectra and the multi-wavelength 
coverage provide us with high quality measurement of several 
galaxy properties. The paper is organized as follows. In Sections 
|2] and [3] we describe the galaxy sample, the galaxy properties 
we are interested in, how the density contrast has been mea- 
sured and all the advantages and drawbacks of this measure- 
ment. In Section|4]we study the redshift evolution of the colour- 
density relation, together with the dependence on environment 
of other galaxy properties (D4000 A break, [OII]/l3727 equiv- 
alent width). In Section |5] we show how stellar mass depends 
on local density, and in Section |6] we disentangle the correla- 
tion among colour, stellar mass and environment. Finally in|2]we 
discuss the scale dependence of environmental effects. Sections 
[H |9] and [10] are devoted to the discussion of our results and to 
the conclusions. Throughout this paper we use the cosmologi- 
cal parameters Q,„ = 0.25, Oa = 0.75, //o = 70 km s"' Mpc"' 
and //() = 100 h. Magnitudes are expressed in the AB system, 
and absolute magnitudes have Hq - 70 incorporated. Moreover, 
"mass" always means "stellar mass", unless differently stated. 

Parallel analyses about galaxy properties and environment 
have been carried out within the sam e data set i n othe r 
works. Using the densi ty fi eld presented inl Kovac et al.l (1201 Oah . 
iBolzonella etal] (l2009l) and lZucca et"ai] (|2P09) study the Galaxy 
Mass Function and Galax y Lum inosity Function as a func- 
tion of environment, Tasc a et aP (|2009) show h ow the lo- 
cal de nsity affects galaxy m orphology and finally Cap uti et al] 
(|2009l) and [S ilverman et alJ (2009a) study the environment sur- 
rounding specific populations (24 jum galaxies and AGN re- 
sp ectively). Moreover, exploiting the group catalog presented 
in iKnobel et al.l (l2009h . [loyinoetal. (2010) and Kovac et al] 



(1201 Obi) study galaxy colou r and morp hology for group and iso- 
lated galaxies. Finallv. .Peng et alj (12010.) draw a global picture 
on the role of mass and environment in influencing galaxy evo- 
lution. 



2. DATA 



2.1. The zCOSMOS "1 Ok-sample" 



zCOSMOS (iLillv et al.ll2007l) is an on-going large spectroscopic 
survey, conceived with the main goal of studying high-redshift 
galaxy environments from very small scales (galaxy groups 
scales) up to larger-scale structures. Two main projects are part 
of the zCOSMOS survey: a bright survey, aimed to measure 
spectroscopic redshifts for ~ 200 00 galaxies with Iab ^ 22.5 
in the ~ 1.7 deg^ COSMOS field dScoviUe et a l. 2007b), and a 
deep survey, with the goal of acquiring ~ 10000 spectra for a 
sample of colour-selected galaxies residing in the redshift range 
1.4 < z < 3.0, in the central 1 deg^ of the COSMOS field. The 
results described in this paper are bas ed on the first ~ 10000 
measured spectra in the bright sample dLillv et alJl2009l) . From 
now on, we will call this sample "lOk-sample". 

Detail s about the zCOSMOS surve y can be found in 
iLillv et alJ (I2007h and iLillv et alj (l2009h In summary spec- 
tra have been collected with flie VIMOS (iLe Fevre et alJl2003h 
multi-spectrograph, at ESO Very Large Telescope. Bright sam- 
ple observations have been performed using the Red Medium- 
Resolution grism (R ~ 600, spectral range 5550-9650 A). The 
final bright sample of ~ 20000 will have an estimated sampling 
rate of ~ 60 - 70% homogeneous over almost the entire field. 



while the 1 Ok-sample has a less homogeneous and lower sam- 
pling rate, on average of ~ 33%, and it actually covers an area 
of- 1.4 deg^. 

Fo r the data reduction we refer the reader to iLillv et all 
(l2009h . Here it is worth noticing that thanks to repeated obser- 
vations carried out on a subsample of galaxies, a measurement 
redshift error of ~ 100 km/s has been estimated. For this paper, 
a subsample of galaxies with an overall 98.4% reliability rate in 
the redshift measurement have been used. This leave us with a 
sample of ~ 7800 galaxies in the redshift range 0.1 < z < 1.0. 

Finally, the COSMOS field has been covered by multi- 
wave length i magin g. Optical and Near-IR data are described in 
ICapak et al] (|2007|) . See references therein for details on obser- 
vati ons in other wavelength s (XMM-Newton, GALEX, VLA) 
and lKoekemoer et all (l2007h for the HST ACS imaging data. 

2.2. Galaxy properties 

In this paper we analyze environmental effects on galaxy prop- 
erties such as absolute magnitude, stellar mass and spectral in- 
dexes. 

Rest frame absolute magnitudes have been computed from 
the best fitting template normalized to each galaxy photome - 
try and redshift, using the code ZEBRA dFeldmann et alj|2006l) . 
The template set is composed by six observed galaxy spec- 
tra (four from IColeman et aP 1980 a nd two fro m Kinne y et al.l 
1996 ). We refer flie reader to O esch et alj (12009) for more de- 
tails. We tested the robustness of our results using al so absolute 
magnitudes derived with the method described in lllbert et al] 
(2005), a nd with t wo different template sets: the library from 
Bruzual & Chariot (2003), hereafter BC03, and stellar popula- 
tion models generated using the PEGA SE2 population synthesis 
code dFioc & Rocca-Volmerangdll997h . We obtained consistent 
results, within error bars, for different choices of template set 
and/or code. 

Galaxy stellar masses have been computed using a SED fit- 
ting technique, with the Hyperz mass code (a modified version 
of Hype rz, see IBolzonella et al] 12000). Different SE D libraries 
have been used: the widely-used BC03 library, the iMarastonI 
(12005, hereafter M05) library and Chariot & Bruzual (2007, 
hereafter CB07, private communication) library. The imposed 
Star Formation Histories are 10 exponentially declining SFH 
with e-folding times from 0. 1 to 30 Gyr, plus a model of constant 
star formation. In this paper we use masses derived with BC03 
libraries, a Chabrier Initial Mass Function (.Chabrier..2003.). and 
dust extinction modeled with the Calzetti's law dCalzetti et al.l 
2000). For further details and comparisons among different li- 
braries see Pozzetti et al. (2009). It is worth noticing that results 
obtained in this paper regarding stellar masses are robust using 
masses derived with different libraries. 

In this work me make use of two spectral features: the 
equivalent width of the [0II]/13727 doublet (£^[011] from now 
on) and the amplitude of the 4000A break in its narrow def- 
inition (£)„4000 from now on), as introduced by iBalogh et alJ 
dl999l) . These are computed by an adapted version of the plat- 
efit code (Lamareille et al. 2 006, see also Tremonti et al. 200J 
and lBrinchmann et al.ll2004J for the original pipeline). This code 
fits the stellar contin uum an d absorpt ion li nes making use of 
the STELIB library dLe Bor gne et alj_ 2003 ) of stellar spectra 
and of the GALAXEV dBruzual & CharlotI 12003 ) stellar popu- 
lation synthesis models. We eliminate from our sample spectral 
measurements affected either by large uncertainties in the wave- 
length calibration, by the superposition with strong night sky 
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Fig. 1. B-band rest frame absolute magnitudes as a function of 
redshift for galaxies of the lOk-sample residing in the central 
area of the zCOSMOS field, in the redshift range considered 
in this paper Two time-evolving luminosity volume limits are 
shown with dashed lines: the red long-dashed line represents 
the fainter limit (Mb - -19.3 - z), for which we are com- 
plete up to z = 0.7, and the orange short-dashed line is the limit 
Mb = -20.5 - z, for which we are complete up to z = 1.0 but for 
which we have enough galaxies only for z > 0.4. 



lines or by stellar continuum subtraction problems. Finally, we 
consider only spectral line measurements with a line significance 
(ratio of the maximum flux of a line to the nns of the continuum 
around it) higher than 1.15. This threshold allows us to exclude 
~85% of fake detections, and still to include the majority of real 
emission lines measurements. 

For D„4000 and [OII]/l3727 related measurements, we con- 
sider for our analysis the redshift range 0.48 < z < 1 .0, as the 
[OII]/l3727 line enters in the spectral range of the adopted grism 
only at z = 0.48 (the D4000A break enters at a relatively smaller 
z, but we want to consider the same z range for the two indica- 
tors). In this redshift range, D„4000 has been reliably computed 
for ~ 4000 galaxies, and we have a successful measurement of 
[OII]/l3727 equivalent width for ~ 2500 galaxies. 

2.3. Density 

The detailed description of the environment par ametrization 
within the zCOSMOS lOk-sample is given in iKovac et aLl 
(1201 Oah . Here we give only a general overview. 

The environment surrounding a given galaxy at comoving 
position r has been characterized by the dimensionless density 
contrast 6(r) = [p(r) - Pm{i')]/Pm(i'), where density is smoothed 
with a given filter. In this formula, Pmir) is the mean density at 
the redshift that corresponds to the comoving position r, while 
p{r) is the local density around the given galaxy. 

We remark that for the environment parametrization not only 
the lOk-sample galaxies has been used, but also the remaining 
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Fig. 2. Comparison of local density computed with volume lim- 
ited tracers ('vol', y axis) and flux limited tracers ('flux', x axis). 
The density computed with the faint and bright volume limited 
tracers is shown in the top and bottom panels, respectively, in 
each case for two redshift bins as indicated on the top of each 
plot. 



~ 30000 galaxies in the full zCOSMOS photometric catalogue, 
exploiting their photometric redshifts. Moreover, boundary cor- 
rections have been applied to density estimates, to take into ac- 
count the fact that for galaxies that reside too close to the field 
edges the filter window will fall out of the survey field. Finally, 
the robustness of this density reconstruction scheme has been 
tested with simulated gala xy catalogue s. 

We refer the reader to ' Kovac et aP (1201 0a') for an overview 
of the properties of the density field reconstructed in the zCOS- 
MOS field. In particular, their figures 22, 23, and 24 show a di- 
rect comparison of the density field an d the distr ibution of galaxy 
groups (Finoguenov et al. 2007; Kno bel et alJ 2009). Another 
analysis of the la rge scale structures in the same field is shown 
in S coville et al.l (l2007a ). For example, they remark that the su- 
permassive structure at z ~ 0.74 seems to be equivalent, in terms 
of mass, to the Coma cluster. 

Given the lower and less homogeneous sampling rate near 
the field boundaries, in this work we limit ourselves to the cen- 
tral area of the zCOSMOS field, enclosed within the rectangular 
area with vertices ra,,,,-,, =149.55, ra^^^ =150.42, dec„i„ =1.75, 
rJpc -2 7 



2.4. Choosing fhe more suitable density computation 
approacii 

The local density p(r) has been computed with a variety of meth- 
ods, including different filters, different tracer galaxy popula- 
tions and different weigh ting schemes for the tracer galaxies. 
While lKovac et al.l ('2010a') discuss all these possibilities, we dis- 
cuss here only the optimal choice adopted for this work. 

The different methods used to estimate the density field pro- 
duce environment parametrizations with significantly different 
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statistical properties (for example the measurement error, or the 
homogeneity across the redshift range explored), and it is impor- 
tant to choose the density estimator best suited for any specific 
scientific analysis. 

For our study, we decided to give priority to the density com- 
puted on as smal l a scale as po ssible. Tests performed using 
mock catalogues (JKovac et al.ll20 10a) have shown that with the 
zCOSMOS lOk-sample we can safely reconstruct the galaxy lo- 
cal environment up to z ~ 1 using cylindrical filters with half 
depth of 1000 km/s, and either a fixed radius R >3 h 'Mpc, or 
a variable radius equal to the projected distance of at least the 
5''' near est neighbou r ('n.n.' from now on). As one can see in 
Fig. 6 of IKovac et al.l (l2010ah . at least at the highest densities the 
distance D to the 5''" n.n. is always smaller than 3 h 'Mpc. Thus 
we will use the density computed with the 5''' n.n. approach. 
However with this method the scale R varies as a function of en- 
vironment, being larger for lower densities. This has to be taken 
into account in the interpretation of our results. 

Moreover, different tracer galaxy populations have been used 
for the density computation: the entire flux limited sample (Iab ^ 
22.5) and two different luminosity-selected volume-limited sub- 
samples, that satisfy the criteria Mb < -19.3 - z and Mg < 
-20.5 - z- The redshift dependence in these equations has been 
chosen to take into account the effects of an average passive evo- 
lution that should include the majority of galaxy types, and the 
intercepts values have been chosen in order to defi ne two com- 
plete s amples at z = 0.7 and z - 1.0, respectively dKovac et al.l 
1201 Oah . This is clearly shown in Fig. \T\ The figure shows the 
B-band absolute magnitude distribution as a function of redshift 
for the lOk-sample, in the redshift range considered in this paper 
The two time-evolving luminosity limits are plotted as dashed 
lines. From the figure, it is evident that the fainter subsample is 
defined to inlcude all the objects at its higher redshift (z = 0.7), 
while the brighter subsample at z = 1 remains about 0.5 mag 
brighter than the faintest galaxies that we can see in the 10k- 
sample. This is due to the fact that our survey is flux-limited in 
I-band, which corresponds to the restframe B-band at z ~ 0.7. 
Thus, our sample is complete in restframe B-band at this red- 
shift. At higher redshift, because of the different evolution of dif- 
ferent galaxy populations, we start missing the early type galax- 
ies at the faintest B-band magnitudes we can reach. The brighter 
magnitude limit that we used ensures us to be complete for all 
galaxy populations. 

The advantage in using the ffux limited tracers is that we can 
exploit our entire data set, and we can measure the density con- 
trast down to the smallest scales adopting the 5''' n.n. method. 
There are two drawbacks in using these tracers: the tracer pop- 
ulation is not uniform with redshift, as it becomes brighter for 
higher z, and the mean distance to the 5''' n.n. increases with z. 
Both these issues have to be considered when interpreting any 
evolutionary environmental effect on galaxy properties. These 
disadvantages are not present if we use a luminosity-selected 
volume-limited sample of tracers. For these reasons, we decided 
to use the density contrast computed with the tracers limited at 
Mb < -19.3 - z or Mb < -20.5 - z, according to the redshift 
range we want to explore. The limit at Mb < -20.5 -z at low z in- 
cludes too few galaxies, so these tracers can be safely used only 
in the range 0.4 < z < 1.0. In contrast, the fainter luminosity- 
limited sample {Mb < -19.3 - z) is complete only up to z = 0.7, 
but can be used starting from z = 0. 1 . In Fig.|2]we show the com- 
parison between the density computed with flux limited tracers 
and volume limited tracers. 

Finally, different weights can be assigned to the tracer galax- 
ies. For example, we can weight galaxies by their stellar mass or 
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Fig. 3. Left panel. Fraction of red (U-B > 1 .0) galaxies as a func- 
tion of the density contrast, in equipopulated density bins, for a 
luminosity-selected volume-limited sample with Mb < -19.3-z 
in the redshift range 0.5 < z < 0.7. Density is computed with 
the 5''' nearest neighbour method, using volume limited tracers 
(Mb < -19.3 - z). The values on the x axis represent the median 
density values in each density bin. Red squares are for density 
based on counts ('CWD', see text), while orange triangles for 
density weighted by stellar mass ('MWD'). Vertical error bars 
are computed applying the binomial formula. Solid line is the 
linear fit of squares, dashed line the linear fit of triangles. The 
slopes of the fits with their errors are quoted in the text. Right 
panel. As in the left panel, but in this case, before computing 
the MWD to each galaxy it has been assigned the entire set of 
properties (mass, colour etc..) of another galaxy near in redshift, 
chosen randomly. This way galaxy positions and redshift are not 
correlated with galaxy properties. See the text for further details. 



their luminosity in a given band, or simply by their number As 
this is a subtle issue for what concern the general definition (and 
thus interpretation) of local density itself, we discuss explicitly 
advantages and drawbacks of different weighting schemes in the 
following Section. 

3. The mass-weighted density and the 
colour-density relation 

It has been shown that galaxy formation is 'biased', in the sense 
that galaxies have been formed earlier in higher matter-density 
peaks (see for example Marinoni et al. 2005 ). This suggests that 
the galaxy total mass, if not even the (dark) matter distribution 
on larger scales, is a key ingredient to understand galaxy forma- 
tion and evolution. Moreover, we know that galaxy stellar mass 
correlates with the dark matter halo mass, for both early an d 
late type galaxies (iMandelbaum et al.ll2006t lYanget alJbOOSi) . 
This leads to the idea that a density estimator that takes into 
account the stellar mass of the tracer galaxies could be more 
physically motivated than simple galaxy counts. From now on, 
we will define these two different environment parametrizations 
as 'mass-weighted density' (MWD) and ' count-weighted den - 
sity' (CWD). We also refer the reader to lKovacetal.1 (l2009l) . 
who study in detail how the reconstructed galaxy density field 
within the zCOSMOS lOk-sample relates to the underlying mat- 
ter density field. 

The concern in using the MWD is the existence of a correla- 
tion between stellar mass and other galaxy properties. The pos- 
sible environmental dependence of these properties risks to be 
spuriously increased by their dependence on mass. A typical ex- 
ample, that moreover is the one we are mostly interested in, is the 
colour-mass relation. The left panel of Fig.[3]shows the fraction 
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of red (U-B> 1.0) galaxies as a function of the density contrast, 
in equipopulated density bins, for the same luminosity-selected 
volume-limited sample used for the computation of the density 
contrast, in the redshift range 0.5 < z < 0.7 (this cut leaves us 
with roughly 1200 galaxies). Density is computed with the 5''' 
nearest neighbour method, using the fainter volume limited trac- 
ers (Mb < -19.3 - z). Orange triangles are for MWD, while red 
squares are for CWD. Vertical error bars are computed applying 
the binomial formula (see Section|4]for details). The dashed line 
is the linear fit of triangles, the solid line the linear fit of squares. 
The slopes of the fits with their eiTors are 0.31 + / - 0.03 and 
0.21 -H / - 0.04 for MWD and CWD respectively. We note that 
the colour-density relation appears to be stronger and more sig- 
nificant using MWD. 

We tested the influence of the underlying colour-mass rela- 
tion on this slope change by shuffling galaxy properties, while 
preserving their coordinates and redshift. The shuffle has been 
performed by assigning to each galaxy the complete set of prop- 
erties (mass, colour, luminosity...) of another galaxy chosen ran- 
domly in the same narrow redshift bin. We are aware that this 
choice can make the randomization not ideal, as galaxy proper- 
ties are not completely un-correlated with galaxy 3D positions, 
but in this way we preserve the redshift distribution of galaxy 
properties. We then recalculated the CWD and MWD using the 
real galaxy positions, but the randomly assigned masses. This 
way, the CWD estimate does not change, but the 3D mass distri- 
bution (as traced by galaxy stellar mass) is not linked any more 
to galaxy position, and thus the MWD estimate is different from 
the original one. 

The right panel of Fig. [3] shows the colour-density relation 
for the same galaxies as in the left panel , but in this case we used 
the new CWD and MWD and the 'shuffled' colours. Symbols are 
the same as in the left panel. As expected, there is no significant 
colour-density relation when using the CWD and the random 
colours, as these colours are completely uncoiTelated with the 
original galaxy position. The measured slope of the linear fit is 
0.06 + / - 0.04. The fit is not totally flat, possibly due to the 
non-optimal randomization. On the contrary, when we use the 
MWD, we do see a colour-density relation (the measured slope 
of the linear fit is 0. 14-I-/-0.03). As we know that the colour is no 
more (or very weakly) related to the galaxy position, but it does 
depend on mass, the difference that we see in the colour-density 
relation when using CWD and MWD is a spurious effect pro- 
duced by the colour-mass relation. Although the mass-weighted 
density may be considered a physically motivated environment 
estimator, this test suggests us that it is safer to use the density 
computed with pure counts. Thus we will use the CWD in this 
work. 

Summarizing our overall choice, in this work we will use a 
density contrast estimator based on the 5''' n.n. method, com- 
puted with volume-limited and 'count-weighted' tracers. 

4. Redshift evolution of environmental effects 

As a first step, we want to exploit the broad redshift range cov- 
ered by the zCOSMOS lOk-sample to investigate the evolution 
with cosmic time of environmental effects on galaxy proper- 
ties. If we could constrain the epoch (if any) at which environ- 
ment begins playing a role in galaxy evolution, and understand 
whether the properties of galaxies residing in different environ- 
ments evolve differently with cosmic time, we could shed more 
light on galaxy evolution mechanisms. 

In this analysis, we want to use the density computed with 
the bright volume limited tracers complete up to z = 1.0 (Mb < 
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Fig. 4. Fraction of red galaxies as a function of the density con- 
trast. Galaxies are defined 'red' when they have (U - B) > 1.0. 
Different symbols are for different density measurement meth- 
ods: red squares for bright volume limited tracers [Mb < -20.5 - 
z), while orange triangles for flux limited tracers. The values on 
the X axis represent the median density values in each density 
bin. The continuous red line is linear fits of the squares, and 
dashed orange line is the linear fit of the triangles. Their slopes 
with associated errors can be found in Table[T] 'Ntot' is the total 
number of galaxies considered in each panel. 



-20.5 -z). As discussed in Section l274l with these tracers we can 
reliably measure the density only for z > 0.4. So we restrict this 
study to the redshift range 0.4-1.0. 

The galaxy set used for this analysis is the luminosity- 
selected volume-limited subsample that corresponds to the 
bright tracers galaxies. As explained in Sect l2.4l the dimming 
of 1 mag per unit redshift interval is an average evolution that 
should include roughly the typical passive evolution of the ma- 
jority of galaxy types. The value of the intercept is chosen to 
have a complete sample at z = 1 . This limit allows us to approxi- 
mately follow an homogeneous and complete galaxy population 
up to z = 1 . With this cut in luminosity, and considering the red- 
shift range 0.4 < z < 1 .0, we are left with a total sample of 11 33 
galaxies (out of 3868 galaxies in the flux limited sample over the 
same redshift range). 

We analyzed the widely studied colour-density relation, fo- 
cusing on the restframe U - B colour. The zCOSMOS bright 
sample is selected in the observed I-band magnitude, that cor- 
responds to a selection in the absolute B-band magnitude at 
z ~ 0.8. This way, our sample is not affected by any signifi- 
cant U-B color incompleteness up the highest redshift we are 
interested in (z ~ 1). 

The plot in Fig. |4] shows the fraction of red galaxies (red 
squares) as a function of the density contrast, in three differ- 
ent redshift bins. Galaxies are defined 'red' when they have 
(U - B) > 1.0. We choose this colour threshold as it is roughly 
in the middle of the so called green valley of the bimodal colour 
distribution. We keep this value fixed at all explored redshifts, as 
we do not see any strong evolution of the green valley location 
with z. The density bins are chosen so that the total sample used 
in this analysis (the considered 1133 galaxies) is divided into 
three equipopulated bins. This way the density bins in each panel 
are not exactly equipopulated, but the environment thresholds 
are constant with redshift. This implies that we are observing, at 
different redshifts, comparable over-/under-densities, disregard- 
ing the average evolution of the density contrast. Vertical eiTor 
bars on the red fraction values are the Icr confidence level for 
the binomial statistics, computed according to the approxima- 
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Fig. 5. Top. Fraction of 'passive' galaxies as a function of the 
density contrast. Galaxies are defined 'passive' when they have 
Z)„4000 > 1.40. Symbols and lines have the same meaning 
as in the plots of Fig. |4] Bottom. Fraction of galaxies with 
log(£'W[(9//]) > 1.15 as a function of the density contrast. 
Symbols and lines as in the top panels. 



tion given in iGehrelsl ( Il986h : cp- — frfhin, where /, is the frac- 
tion of 'red' galaxies, ft is the fraction of 'blue' ones (= 1 - /,) 
and n is the total number of galaxies in the density bin. When 
the fraction for which we want to compute the Icr confidence 
level is zero (/one), we fix its lower (/upper) limit to zero (/one), 
and its upper (/lower) Icr level at o" = 0.5 /n, as suggested by 
iDe Propris etal] ( |2004|) . 

We fit the colour-density relation with a linear fit (red solid 
line), taking into account the vertical error bars of the points. We 
disregarded the estimated horizontal error on the median density 
value, as it is so small (on average around 0.06 dex) that the error 
on the linear fit is dominated by the error on the red galaxies 
fraction. In contrast, we verified that changing the number of 
density bins (compatibly with the total number of galaxies) does 
not change significantly the result of the fit. These considerations 
hold for all our subsequent analyses. 

The derived slope of the linear fit and its error are quoted in 
the second column of Table[T] The significance associated to the 
colour-density trend is ~ 3cr, ~ 3cr and Icr from the lowest to the 
highest redshift bin. The colour-density relation, present at low 
redshift, is still present up to z = 1 . Moreover, there is a trend of 
weakening of the colour-density relation for increasing redshift, 
as the slope of the fit becomes less steep at higher redshift. 

As a reference, we list in Table [T] the slopes (and their er- 
ror bars) of the linear fit of the colour-density relation that we 
find when using also different tracers of the local density (vol- 
ume limited and MWD, flux limited and CWD, and flux limited 
and MWD, in the third, fourth and fifth column, respectively). 
The reader can notice the foreseen effect of a enhanced colour- 



density relation when 1) investigating the smaller scales probed 
by flux limited tracers and 2) using the mass-weighted density. 
The effect given by the flux limited tracers sample is evident at 
least for 0.4 < z < 0.7, where the flux limited sample is much 
fainter than the volume limited one (see Fig.[T]). 

Usually, galaxy colours can be considered as a proxy for star 
formation activity, with redder galaxies having less on-going 
star formation. We thus analyze the dependence on environ- 
ment of two other star formation indicators: the D„4000 and 
the ^^[OII]. In particular, D„4000 is quite well related with 
the U - B colour, as the U and B bands bracket the break it- 
self. Nevertheless, the D„4000 is a direct measure of the stellar 
population, and it is less aff'ected by reddening or uncertainties 
due to the SE P fitting necessary for the colour computation (see 
for example in lFranzetti et alj|2007h . 

We therefore identify two logical equivalents of our 'red' 
galaxy sample by defining a 'passive' galaxy subsample, consist- 
ing of all galaxies having Z)„4000 > 1.4, and an 'active' galaxy 
subsample including th ose galaxies with \o ^{E W\OIIYi > 1.15. 
The reader is referred to lMaier et al.l (l2009h and lSilverman et al] 
(2009b) for more details on the star formation activity within 
zCOSMOS galaxies. 

These thresholds correspond roughly to the median values 
of the total D„4000 and £^[011] distributions. As for t/ - B 
colour, the D„4000 threshold does not evolve significantly with 
redshift, thus we used it for all the redshift bins investigated; on 
the contrary, the £^[011] mean value increases with redshift, but 
we verified that the dependence of the 'active' galaxy fraction 
on environment does not change significantly if we substitute 
the fixed threshold with an evolving one (with the exception of a 
different normalization). 

In Fig. |5] we show the fraction of 'passive' (top row) and 
'active' galaxies (bottom row) as a function of density, for three 
different redshift bins (the three columns). The symbols and the 
solid lines in these panels have the same meaning as those in Fig. 
m Note that these features enter the observed wavelength range 
at about z - 0.48, thus the lower limit of the first redshift bin 
is different from the one in Fig. |4] The slopes of the linear fits, 
together with their Icr levels, are quoted in Table [T] 

Figure |5] shows that the environmental effects on D„4000 as 
a function of redshift mirror those on galaxy U - B colour, with 
comparable slopes and normalizations (inside the error bars). As 
the 'red' galaxies, the 'passive' ones reside preferentially in high 
density regions, and this D„4000-density relation is visible in the 
entire range 0.48 < z < 1, weakening for increasing redshift. On 
the contrary, the fraction of 'active' galaxies is higher in low 
densities. Even if less significant than the detection of D„4000- 
density relation, we find that the fiW [Oil] -density relation holds 
up to z ~ 1, and like the others it becomes weaker for increas- 
ing redshift. In Table [T] we quote the slopes of the linear fits of 
these relations also when MWD and flux limited tracers are used 
to estimate the environment. As for the colour-density relation, 
results are compatible within errors irrespectively of the density 
estimator used, with a general indication for stronger and more 
significant trends in the MWD case. This is mainly due to the 
relation between stellar mass and both D„4000 and £^[011], as 
discussed in Section [3] Moreover, it seems that above z ~ 0.7 
the £^[011] is less sensitive to environment than colour and 
Z)„4000. A more detailed analysis of this issue is deferred to 
future work. 

We remark that we excluded from this analysis all galaxies 
with low quality spectral features measurements as described in 
Section l272l These galaxies are ~ 15% and ~ 45% of the sample 
used for the colour-density analysis, considering the D„4000 and 
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Table 1. Slopes and their Icr confidence levels of the linear fits 
shown in Fig.|4]and|5](second column), that is the slopes of the 
'property' -density relations when density is computed with vol- 
ume limited tracers (Mb < -20.5 -z) and 'count-weighted' den- 
sity (CWD), with 'property' standing for U - B colour, D„4000 
and ^^[OII] as indicated in the first column. From the third to 
the fifth column the same slopes are quoted, but in these cases 
they are obtained with density computed with different tracers 
(volume limited and MWD, flux limited and CWD, flux limited 
and MWD, respectively). 



z range 


volume limited tracers 


flux limited tracers 




CWD 


MWD 


CWD 


MWD 


U-B 










0.40 - 0.7 


0.23+0.08 


0.33+0.07 


0.26+0.07 


0.44+0.06 


0.7 - 0.85 


0.23+0.08 


0.25+0.06 


0.18+0.08 


0.25+0.06 


0.85 - 1.0 


0.14+0.07 


0.19+0.06 


0.15+0.07 


0.20+0.06 


Z)„4000 










0.48 - 0.7 


0.30+0.09 


0.38+0.08 


0.38+0.08 


0.48+0.07 


0.7 - 0.85 


0.22+0.08 


0.23+0.07 


0.19+0.08 


0.28+0.06 


0.85 - 1.0 


0.13+0.08 


0.15+0.07 


0.11+0.07 


0.16+0.06 


EWIOII] 










0.48 - 0.7 


-0.20+0.11 


-0.27+0.09 


-0.28+0.09 


-0.35+0.09 


0.7 - 0.85 


-0.07+0.11 


-0.05+0.08 


-0.06+0.10 


-0.13+0.08 


0.85 - 1.0 


-0.06+0.09 


-0.11+0.09 


-0.06+0.09 


-0.10+0.08 



the £^[011] measurement respectively. To be sure that this spec- 
tral quality selection does not introduce any bias as a function of 
density, we re-computed the D„4000- and £^[011]- density re- 
lations considering also these previously excluded galaxies. We 
found results compatible (within errors) with the ones shown in 

Fig.ia 

Summarizing, in Fig. |4] and |5] we have shown how differ- 
ent galaxy properties depend on environment within a complete 
luminosity-selected volume-limited sample, and how these de- 
pendences evolve with redshift. Globally, these results confirm 
already known trends: the fraction of red galaxies is higher in 
high densities, and the colour-density relation becomes weaker 
for increasing redshift. We have also shown that the fraction of 
galaxies with high D„4000 depends on environment in a way 
that is very similar to the U-B colour (as both red colour and 
high D„4000 are indicators of old and passively evolving galax- 
ies), and this is the first time that this is shown up to z ~ 1. 
Conversely we found that the 'active' galaxies (those with larger 
£^^[011]) reside preferentially in low densities. As for colour 
and D„4000, the £'W[OII] -density relation holds up to z ~ 1, al- 
though with less significance, and weakens for increasing red- 
shift, further confirming the general trends observed with the 
other star formation indicators. We refer the reader to Section 
[8]for the comparison of these findings with previous works, and 
to Section|9]for the discussion of their implications. 

The comparative study of the environmental dependences of 
different galaxy properties at different epochs is of crucial im- 
portance to shed more light on galaxy evolution. We defer this 
analysis to future work, while here from now on we will explore 
more in details only the colour-density relation. 



5. Mass segregation as a function of environment 

It has akeady been shown that galaxy stellar mass (M) does 
depend on environment, at least weakly, at both z ~ 0.1 
and up to z = 1 .4 (Kauffmannet al. 2004; B undv et al.ll2006t 
IScodeggio et al.l2009t:iBolzonella et al.i2009i) . There is still con- 
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Fig. 6. The mass-redshift plane for galaxies of the lOk-sample 
residing in the central area of the zCOSMOS field, with red 
curve representing the 95% mass completeness described in 
Section |5] The three redshift bins considered in the analysis and 
their chosen mass limits are indicated with straight lines. 



siderable debate whether this mass-density relation is the driver 
of environmental effects upon other properties of a galaxy such 
as its colour and efficiency of star formation, given their known 
correlations with stellar mass. For example, in the local uni- 
verse galaxy colour, D„4000 and sSFR are found to depend on 
environment even when mass is fixed dKauffmann et al.l 120041 : 
Baldrv et al. 2006), but this has not been found in th e rang e 
0.2 < z < 1.4 within VVDS data (IScodeggio et all 120091) . 
Although Scodeg gio et al.l (12009.) themselves explain the possi- 
ble reasons for this discrepancy, up to now no other study of this 
kind has been performed at intermediate-high redshift. With the 
zCOSMOS lOk-sample data we can now add a new insight on 
this issue in the range 0.1 < z < 1.0. 

The detailed analysis of the mass-d ensity relation wi t hin th e 
zCOSMOS data is fufly described in iBolzonella eTaP (l2009l) . 
where the Galaxy Stellar Mass Function (GSMF) is studied for 
different galaxy subsamples and environments. For the sake of 
clarity and completeness, in this Section we will simply present 
how we select a mass limited subsample suited for our follow- 
ing analysis, and how the stellar mass depends on environment 
in particular within this selected sample. Then, in Section |6] we 
will disentangle the triple relation among colour, stellar mass 
and local density. 

As in all flux limited samples, the range of luminosities 
and masses allowed in our sample varies with redshift (see Fig. 
[U and |6]l. Given the scatter in the luminosity-mass relation, a 
luminosity-selected volume-limited sample at a given redshift 
will not be complete in stellar mass. As we now require a mass- 
complete sample for our analysis, we can not use the luminosity- 
selected volume-limited sample described in SectionH) Since the 
mass-to-light ratio (M/L) is different for different galaxy types, 
to select a mass-complete sample we need to identify at any red- 
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shift the smallest observable stellar mass for the red/early type 
galaxies, which have on the mean the higher M/L at any given 
mass. In practice, as proposed by Zucca et al. (2006), we fit the 
Spectral Energy Distri bution of our galaxie s to six templates 
(four observed spectra, Coleman etal.|[l980l and two starburst 
SEDs. iBruzual A. & Charloall993h . and we compute our mass 
completeness using those galaxies that have been assigned the 
earliest template, i.e. the E/SO template. We define as the limit 
mass the mass at which we are 95% complete for these galaxies, 
at any given redshift. Figure |6] shows the mass-redshift plane of 
the lOk-sample, the red curve representing the 95% mass com- 
pleteness described above. 

As it can be noticed, the mass limit at z = 1 .0 is very high. If 
we used it in the whole range 0.1 < z < 1 .0 in order to select a 
mass limited sample homogeneous at all redshifts, we would be 
left with too few galaxies at low z- We will also show that for the 
analysis carried on in this part of our work we are interested in 
as large a mass range as possible at any redshift. For these two 
reasons it is more convenient to separately adopt the smallest 
limit mass we can reach in any redshift bin. Within each z bin, 
we keep the limit mass constant (unlike in the case of the limiting 
magnitude), as it has been shown that the Galaxy Stellar Mass 
Function does not ev olve significantly in the redshift range we 
are considering (e.g. . iPozzetti et al.ll2007h . 

To cover the largest available mass range we start our anal- 
ysis from z = 0. 1 . This choice prevents us from using in the en- 
tire redshift range the density contrast computed with the bright 
volume limited tracers, as it is reliable only for z > 0.4. In con- 
trast, the fainter volume limited tracers (Mb < -19.3 - z) rep- 
resent a complete sample only up to z - 0.7. We decided to 
use the fainter volume limited tracers in the redshift interval 0. 1- 
0.7, and the bright one for the remaining bin 0.7 < z < 1.0. 
We verified that this choice does not introduce any spurious de- 
pendence on redshift, because of the change, at z = 0.7, of the 
tracers used to measure the density. We took as a reference the 
redshift bin 0.5 < z < 0.7, for which we have the density com- 
puted with both the samples, and we performed all the analyses 
shown in Sections |5] and |6] with both the tracers. We find that 
the results are always consistent using the two tracers samples, 
as shown for example in Fig. Q and [8] (see the text below for 
th e details about these fig ures). The same agreement is found 
in lBolzonella et all (l2009l) . where the galaxy stellar mass func- 
tions in low and high density environment are computed for both 
tracer types. Therefore we believe to be safe the comparison of 
our analysis at z < 0.7 and z > 0.7 using fainter and brighter 
tracers, respectively. 

In summary, we selected the three redshift bins [0.1- 
0.5], [0.5-0.7], [0.7-1.0], with the mass Hmits Mu,„ given by 
log(M/,v„/M0) = 10.2, 10.65 and 11.05 respectively. We then 
studied the mass-density relation within these mass limited sam- 
ples. This is shown in Fig. |7j where we plot the median stel- 
lar mass (red squares and blue diamonds) as a function of den- 
sity. Red squares are computed using the fainter volume lim- 
ited tracers, blue diamonds using the bright ones. The density 
bins are chosen in the following way. For the fainter tracers, 
we selected all the galaxies in the range 0.1 < z < 0.7 with 
log(M/MQ) < 10.65 (our highest mass limit for z < 0.7), and 
then we divide this sample in 4 equipopulated density bins. Then 
we use the so-found density thresholds in the two redshift bins 
0.1-0.5 and 0.5-0.7. Similarly, for the brighter tracers we took 
all the galaxies with log(M/M0) < 11.05 within 0.5 < z < 1.0, 
we divide this sample in 4 equipopulated density bins, and we 
use these density thresholds in the two redshift bins 0.5-0.7 and 
0.7-1.0. The red lines in Fig. [7] are linear fits to the squares, the 
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Fig. 7. Median stellar mass (red squares and blue diamonds) as 
a function of the 'count-weighted' density contrast with volume 
limited tracers, for the three mass limited subsamples described 
in Section |5] Black points are single galaxies. Red squares 
are computed using the fainter volume limited tracers (Mb < 
-19.3-z), blue diamonds using the bright ones (Mb < -20.5-z). 
X axis values for squares and diamonds represent the median 
density values in each density bin. Vertical error bars associated 
to median stellar mass values are obtained with bootstrapping 
technique. Black thick(/thin) lines are the 15% and the 85% of 
the mass distribution in each density bin, for the faint(/bright) 
volume-limited tracers. Horizontal dashed lines show the mass 
limits in each redshift bin. The red straight lines are linear fits 
of the squares, the dashed blue lines of the diamonds, the slopes 
and their errors being quoted in the text. The total number of 
galaxies considered in each redshift bin is quoted in the upper 
part of each panel. 



blue dashed lines to the diamonds. The slope values and their 
errors are 0.9 + 0.03 and 0.8 + 0.04 for the squares in the first 
and second redshift bin, and 0.8 + 0.04 and 0.6 + 0.03 for the 
diamonds in the second and third redshift bin, respectively. 

It can be noticed that the median mass depends on environ- 
ment at z < 0.5, but this dependence weakens for higher red- 
shift. This may not be a pure time-evolution effect, because also 
the mass limit changes with redshift. Considering only the mass 
limited sample in 0.1 < z < 0.5, we verified that within this sam- 
ple the mass-density relation weakens when increasing the mass 
limit. These findings are in agreement with Bolzonel la et ap 
(120091) . who show that the GSMFs in the lOk-sample in the low 
and high density regions differ mainly in the low-mass part, with 
this difference decreasing for higher redshift. 



6. Are the environmental effects on galaxy colour 
driven by mass? 

In the previous Sections we have shown how galaxy U - B 
colour and stellar mass depend on environment. It is also well 
known that colour depends on stellar mass, thus we want to study 
whether the environmental dependence of U-B is simply the ef- 
fect of the combined colour-mass and mass-density relations. We 
address this issue studying the threefold relation among density 
contrast, stellar mass and colour In the following subsections 
we study the colour-density relation in the mass limited samples 
(Section l6Tt and then we discuss how to disentangle the colour- 
mass-environment relation (Section|i 
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Fig. 8. Fraction of red (U - B > 1.0) galaxies as a function of 
density for three mass limited galaxy samples (mass lower limit 
in the label). Red squares are computed using the fainter vol- 
ume limited tracers {Mb < -19.3 - z), blue diamonds using the 
bright ones (Mb < -20.5 -z). The values on the x axis represent 
the median density values in each density bin. The red straight 
lines are linear fits of the squares, the dashed blue lines of the 
diamonds, the slopes and their errors being quoted in the text. 
The total number of galaxies considered in each redshift bin is 
quoted in the lower part of each panel. 



6.1. The colour-density relation in mass-limited subsamples 

The three panels of Fig . [8] show the fraction of red (U - B > 1.0) 
galaxies as a function of density for the three mass limited sam- 
ples presented in Section |5] Red squares are for fainter volume 
limited tracers, and blue diamonds for the brighter ones. The 
density bins along the x axis are selected as in Fig.|7] The linear 
fits to the squares and diamonds are represented with solid and 
dashed lines, respectively. For these fits, we find the following 
slopes (and 1 cr confidence level): 0.17 + 0.04 and 0.05 + 0.05 for 
the squares in the first and second redshift bin, and 0.05 + 0.06 
and -0.01+0.05 for the diamonds in the second and third redshift 
bin, respectively. We find that there is a clear colour-density re- 
lation within the mass-limited sample in the lowest redshift bin, 
but we do not detect it at z > 0.5. We remind the reader that 
this different behaviour at different redshift may not be a pure 
evolutionary effect, because in this figure we consider samples 
with different lower mass limits (see also the last paragraph of 
Sect. |5]). We do not want to make any statement about galaxy 
evolution inspecting Fig. [8] 

It can be noticed that results presented here are quite dif- 
ferent from the ones in Fig. H) where a luminosity-selected 
volume-limited sample was used and where the colour-density 
relation was steeper and with a lower normalization. To under- 
stand these differences we inspect Fig.|9] For the three z bins of 
Fig. [8] it shows the colour-mass plane for all the galaxies (or- 
ange points), with green points representing the galaxies consid- 
ered in the luminosity-selected volume-limited sample of Fig.|4] 
(Mb < -20.5 - z) and the vertical lines the mass limits of Fig. [8] 

The higher normalization in Fig. [8] is explained by the fact 
that with the chosen mass limits we select an higher fraction 
of red galaxies. This effect increases with z, up to the highest 
redshift bin, where almost only red galaxies are selected. This 
observation suggests that the fixed colour cut at U - B = 1 to 
define red galaxies is not well suited for the study of our mass- 
selected samples, and we will expand on this point in the next 
Subsection. 



We also notice how the colour-density relation for the mass- 
selected sample within 0.5 < z < 0.7 (middle panel in Fig. [8]) 
is significantly flatter than the one obtained for the luminosity- 
selected volume-limited sample within roughly the same red- 
shift interval (0.4 < z < 0.7, first panel of Fig. |4|. Observing 
the second panel in Fig. |9l we see that the explanation for 
this flattening is related to the different galaxy populations that 
are selected using luminosity or mass limits. In the luminosity- 
selected volume-limited sample there is a population of lower- 
mass blue galaxies that is not present in the mass limited sam- 
ple due to their low M/L. Moreover, this population has no red 
counterpart (at the same mass) in the luminosity selected sample, 
as these redder galaxies are too faint to satisfy the luminosity- 
selection criterion (they have on the mean an higher M/L). 

This tail of less massive blue galaxies that we find in the 
luminosity-selected sample inhabits preferentially low density 
regions. This can be deduced noticing that the red fractions in 
Fig. |4] are on average always lower than those in Fig. [8] but the 
percentage difference is much larger in the low density bins. It 
remains to be clarified whether these galaxies reside mainly in 
low density environment only because of their low mass (the 
mass segreg ation effect that we showed in Fig. |7] and that is 
detailed in .BolzoneUa et alll2009l) . or because environment also 
plays a role in shaping their colour. We will disentangle this de- 
pendence in the following Subsection. 

Here, it is worth remarking the following. We are aware that 
the M/L spread is particularly large in our sample because of 
the relatively blue band-passes we use (U- and B-band). These 
bands are more sensitive to dust reddening, and moreover the 
light at these wavelengths is dominated by younger stars, that 
not only are a small fraction of the mass, but also can be the 
product of transient phenomena, such as sudden bursts of star 
formation. As already explained, our choice of the B-band for 
the luminosity selection and of the U-B colour for the analysis 
is dictated by the fact that the 1-band flux limit of our survey 
(Iab ^ 22.5) corresponds at the rest frame B-band at z ~ 0.8. 
This allows our sample to be not significantly affected by colour 
incompleteness in the redshift range we are interested in (z < 1). 
Moreover, in a luminosity-selected sample the M/L spread is 
in general unavoidable, irrespectively of the band used for the 
selection, because it is produced at least by the variety of star 
formation histories. This gives rise in turn to a large spread in 
stellar masses, and together with the colour-mass relation, this 
causes the colour-density relation in a luminosity-selected sam- 
ple to be strongly biased by the mass-density relation. In a mass- 
limited sample this bias is weaker, as the allowed mass range is 
smaller, but the mass-density relation does not disappear. Thus, 
if we want to study the direct environmental effects on galaxy 
properties, such as the colour, without being biased by the mass 
segregation as a function of environment, we need to select our 
sample(s) in narrow mass bins. 

6.2. Disentangling the colour-mass-density relation 

Following the discussion detailed above, in this Section we study 
the colour-density relation in mass bins, in order to understand 
whether it is only the result of combining the mass-density and 
colour-mass relations. First, we divide the three mass limited 
samples in mass bins of A log(M/M0) = 0.2, in order to can- 
cel the mass-density relation. We verified that in such narrow 
mass bins the stellar mass does not depend any more on en- 
vironment. Secondly, as previously discussed, the colour cut at 
U - B - 1 is not well suited to define 'red' galaxies in a sample 
with log(M/M0) > 10.7, thus we introduce a new definition of 
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Fig. 10. Colour-density relation using number density and volume limited tracers for subsamples of galaxies divided in three redshift 
bins (columns, redshift bins quoted on the top) and five mass bins (rows, mass bins quoted on the right, in logarithmic scale). For 
z < 0.7 we use the fainter volume limited tracers {Mb < -19.3 - z), for z > 0.7 brighter ones (Mb < -20.5 - z). Red circles show 
the fraction of galaxies with U - B colour redder than a threshold that depends on the stellar mass, as shown by the diagonal red 
line in the first panel of Fig.[TT](see also text for details). The values on the x axis are the median density values in each density bin. 
Red lines represent the linear fit of circles. Their slopes and associated Icr confidence levels are quoted in Table |2] In each panel, 
'Ntot' is the total number of galaxies in that mass and redshift bin, while the number of galaxies satisfying the colour threshold is 
quoted in Table |2] 

Table 2. Slopes of the linear fits plotted in the panels of Fig. [TO] Different rows are for the different mass ranges, and the columns 
for the three redshift bins. A^,, ,. is the number of 'very red' galaxies, and Nu,/ the total number of galaxies in each panel. 



Mass range 



0.1<z<0.5 



0.5< z <0.7 



0.7<z<1.0 



Slope±lo- N,,r/N,o, Slope±lcr N,„JN,„ Slope±lcr A^,,,,/A'„ 



0.06±0.11 42/116 
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Fig. 9. The colour-mass plane for 3 different redshift bins. 
Orange points are all the galaxies in the redshift bin considered, 
green points are those galaxies that respect the luminosity limit 
as in Fig. |4] and vertical lines represent the mass limit at each 
redshift. 



'red' galaxies. We take a colour threshold roughly parallel to the 
red sequence in the colour-mass plane, and ~ 0.3 mag bluer than 
the red sequence itself. We used as reference the red sequence in 
the first redshift bin of Fig. |9] This colour threshold depends on 
the stellar mass, as the red sequence becomes redder for higher 
masses. The first panel of Fig. [TT] shows the colour-mass plane 
within 0.1 < z < 0.5. The red line is the chosen colour threshold. 
As it is always redder than U - B - 1, we define as 'very red 
galaxies' those galaxies with a colour redder than this threshold. 

In Fig. [To] we plot the fraction of 'very red' galaxies (red cir- 
cles) as a function of the density contrast. The columns are for 
three redshift bins as indicated on top (the same as in Fig. |8]l, 
while the rows represent different mass bins: from top to bottom 
we consider the bins 10.25 - 10.45, 10.45 - 10.65, 10.65 - 10.85, 
10.85 - 1 1 .05 and 1 1 .05 - 1 1 .25 in log{M/Mo) units, as indicated 
in the labels on the right of the right-most panels. The lowest 
mass considered in each redshift bin is always equal or bigger 
than the mass limit at that given z- For this analysis we used 
CWD and fainter volume limited tracers for z < 0.7, and CWD 
and brighter volume limited tracers for z > 0.7. We verified that 
our results in the range 0.5 < z < 0.7 do not change using fainter 
or brighter volume limited tracers. The dashed-dotted line is the 
linear fit of the circles (considering their error), and 'Ntot' indi- 
cates the total number of galaxies in each panel. Table |2] shows 
the fit slopes with their error, together with the number of 'very 
red' galaxies. 

Figure[TO]shows that the fraction of 'very red' galaxies seems 
generally not to depend on environment once mass is fixed, but 
we can notice some exceptions. In fact, in the redshift range 
0.1 < z < 0.5 we see that the fraction of 'very red' galaxies 
depends on environment for log(M) < 10.7. This fraction seems 
to show some trend with density also for the highest masses ex- 
plored (log(M/M0) > 1 1 .0), but for these high masses we have 
larger error bars due t o the lower s t atistics . We note that the same 
trend is observed by iTasca et al] (|2009|) . when they study the 
fraction of the more massive early type galaxies as a function 
of local density. We will explore this mass regime with higher 
statistics with the complete zCOSMOS bright sample. 

We addressed this issue more in detail, focusing our attention 
on the redshift range 0.1 < z < 0.5, where the colour-mass- 
density relation can be better studied thanks to the larger mass 
range spanned above the mass limit. The aim is to analyze in a 
more continuous and smooth way the colour-density relation as 
a function of mass bins. 



The first panel of Fig. [TT] shows a zoom-in of the first panel 
of Fig. [9] but in this case only above the completeness mass limit 
log(M/M0) = 10.2. The black solid curves are the density con- 
tours and the green dashed lines are the 15''', 50''' and 85''' per- 
centiles of the colour distribution in not-independent mass bin 
of Alog(M/MQ) - 0.2, each shifted from the previous one by 
A.5/,,// log(M/M0) = 0.05. In this panel we show with a red di- 
agonal line the colour threshold used also in Fig. [10] to define 
'very red' galaxies. For each mass bin, we computed the frac- 
tion of 'very red' galaxies in four density bins, defined as in Fig. 
[To] In the second panel of Fig. [TT]we plotted the very red galaxy 
fraction for the lowest and highest density bins, with squares and 
circles respectively. Then in each mass bin we computed a linear 
fit of the four fraction values as a function of density, and in the 
third panel of Fig. [TT] we plotted the significance (the cr level) 
of the slope of these fits. The horizontal dashed line is the 3cr 
significance level, for reference. 

Figure [TT] more clearly shows the results of Fig. [10] If we 
fix the galaxy stellar mass, some colour-density relation is still 
measurable at least at a 3cr significance level for galaxies with 
log(M/M0) < 10.6. We note that for these low masses the galaxy 
population is still composed by two main colour sequences, or 
clouds, divided by the so called green valley. On the contrary, ac- 
cording to the colour-mass relation, galaxies with higher masses 
have mainly redder colours. 

The central panel of Fig. [TT] indirectly shows also that a 
colour-mass relation holds irrespectively of environment. In fact 
the definition of 'very red' galaxies implies that we select redder 
and redder galaxies for higher masses. We see that the fraction 
of 'very red' galaxies in high densities is almost constant as a 
function of mass and this means that in this environment galaxy 
colour becomes globally redder with increasing mass. The red- 
dening of the galaxy population for higher masses is even more 
evident in low density regions, where the fraction of 'very red' 
galaxies increases as a function of stellar mass. 

Finally, we observe that the fraction of 'very red' galaxies 
seems to have some (weaker) dependence on the density contrast 
also for the highest masses explored. Unfortunately the number 
of such massive objects in our sample is quite low (less than 100 
galaxies in each of the last three mass bins of Fig. fTTT i. and so 
the observed relation between density and the 'very red' galaxy 
fraction is not significant (if real). It will deserve a deeper in- 
vestigation, once the complete zCOSMOS bright sample will be 
available (almost doubling the statistics of the lOk-sample). 

7. Scale dependence of environmental effects 

In the previous Sections we have studied environment on the 
shortest scales on which we could compute it, given the limits 
imposed by the mean interparticle separation of the galaxies in 
our sample, and thus the limits of density reconstruction reliabil- 
ity found using mock catalogues (Kovac et al. 2010a). Moreover, 
we used volume limited tracers for the reasons listed in Sect. 12.41 
while flux limited t racers can probe shorter scales (see Fig. 13 in 
iKovac et al.ll2010ah . 

The issue of the scale on which density is computed is yet 
of great interest. For example, a still open question is whether 
the eff'ects seen on large scales are only 'residual' of those 
seen at smaller scales, or whether they add information. Some 
works at z ~ 0.1 carried out wifli SPSS and 2dFGRS data 
(iKauffmann et al.l2004tlBlanton et al.l2006l) assess that environ- 
mental dependence of galaxy properties on scales larger than 1 
/z^'Mpc are driven only by the influence of scales < 1/z ' Mpc. 
It could be interesting to extend this study at higher redshift. 
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Fig. 11. Upper panel: colour-mass plane for galaxies in the redshift range 0.1 < z < 0.5, above the mass limit logCM/Mo) = 10.2. 
Green dashed lines represent 15%,50% and 85% of the color distribution, in bins of 0.2 in mass (logarithmic scale, as represented by 
the horizontal line in the bottom right corner). Black contours are the isodensity contours of the plotted points. The red straight line 
is the colour thresholds used to define 'very red' galaxies. It is a straight cut following roughly the red sequence and 0.3 magnitudes 
bluer than the red sequence itself. Middle panel: in not-independent bins of 0.2 in log(M/M0) units, thick circles are the fraction 
of 'very red' galaxies in the highest density quartile and thin squares in the lowest density quartile. The horizontal dashed line at 
fraction=0.5 is for reference. Bottom panel: the significance of the slope of the linear fit of the red fraction as a function of density, 
in each mass bin. In these panels, density is computed with fainter volume limited tracers and without mass weights. 



to provide clues for any interpretation of their results, and also 
possibly to add useful information for planning future survey 
strategies. 

W e thus tried to add ress this issu e using our av ailable data. 
While iKaufFmann et alJ (I2004.) and Blanton et al.l (i2006 ) used 
cylindrical filters with fixed radii, within zCOSMOS lOk-sample 
we can not reach scales as small as ~ l/i"' Mpc with a fixed aper- 
ture, as explained in Section l23] This is expected due to the high 
mean interparticle separation of the lOk-sample, that increases 
from ~ 3 /i^'Mpc to ~ 7.5 /i^'Mpc in the range 0.1 < z < 1.0. 
We do reach scales smaller than l/i ' Mpc only with the S''' n.n. 
technique, using flux limited tracers, and only when considering 
very high densities. 

To b e more precise, as one can see in Fig. 13 in lKovac et al.l 
(l2010ah . using flux limited tracers and the 5''' n.n. method we 
reach scales below l/i"' Mpc for overdensities (log(l + 6) >Qi) 
only within z < 0.5. 

Of course using an adaptive scaling (5''' n.n. ) as a refer- 
ence is not trivial, as scale varies with density. Nevertheless, 
we addressed this issue following an approach s imilar to the 
one shown in Fig. 14 of iKaufimann et all (l2004l) . This figure 



shows that the median value of £)„4000 depends on the number 
of galaxy neighbours, on both small (l/i^'Mpc) and large (5/i"' 
Mpc) scale. The question is whether the density on the larger 
scale still has effects on galaxy D„4000 once the density on the 
smaller scale is fixed. Kauffmann et al. (2004) show that median 
Z)„4000 does not depend on large scale environment, once the 
small scale density is chosen to be enclosed in a small range of 
very low densities or very high densities. We reproduce this anal- 
ysis within the zCOSMOS lOk-sample. We consider the density 
computed on a 5h^^ Mpc Gaussian filter as the large scale envi- 
ronment, and we use the 5 ''' n.n. as the density on the smaller 
scale, of course paying attention at the variation of this scale 
itself with density. We will call these densities [5g] and [5th] re- 
spectively. Our result is shown in Fig. [121 The Figure is divided 
in 4 panels. 

The top left panel shows the fraction of red galaxies (U-B 
> 1.0) for the subsample of galaxies within 0.1 < z < 0.5 and 
with Mb < -18.5 - z, as a function of the density contrast. The 
density is computed with flux limited tracers, and without mass- 
weights, for both the [5th] (red diamonds) and the [5g] (blue 
triangles). The value of the slope of the colour-density relation 
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Fig. 12. Top-left panel. Colour-density relation for galaxies with Mb < -18.5 - z and 0.1 < z < 0.5. We use CWD computed with 
flux limited tracers. Red diamonds are for the 5''' n.n. density estimator ([5th]), while the blue triangles for the density computed 
with a 5 Mpc Gaussian filter ([5g]). x axis values are the median value in four roughly equipopulated bins of the density distribution, 
vertical error bars are given by the usual binomial formula. 'Ntot' is the number of galaxies in the subsample, while the linear fit 
slope value is quoted with its error for both the densities. Top-right panel. Scatter plot of [5th] (x axis) and [5g] (y axis) densities. 
The straight line is the bisector. Bottom-left panel. For any small bin of [5th] density {x axis, equipopulated bins of 300 galaxies), we 
plot the fraction of red galaxies in the first and last quartile of the [5g] density distribution (thick symbols are for the highest density 
bin, light symbols for the lowest density bin). Bottom- right panel. The opposite of the bottom-left panel. For any small bin of [5g] 
density (x axis), we plot the fraction of red galaxies in the first and last quartile of the [5th] density distribution (same symbols as in 
the bottom-left panel). 



is compatible within the two density estimators (scales), but with 
the [5th] we reach higher densities, where we observe the highest 
red fraction. 

The top right panel shows the relation between [5th] and [5g] 
densities, on the x and y axis respectively. The solid line repre- 
sents the bisector, for reference. It can be noticed that for the 
lowest densities, the two environment estimators are very simi- 
lar. This is due to the fact that the distances from the 5''' n.n. are 
getting closer to the scale of the [5g]. On the contrary, for high 
densities, the [5th] density contrast reaches values well above the 
maximum spanned by [5g], because of the smaller and smaller 
scales probed by the [5th]. 

In the bottom panels we want to test whether some 
environment-colour relation is still present on larger scales 
([5g]), when selecting naiTow bins of densities on small scales 
([5th]), and vice-versa. 

The bottom-left panel of Fig. [12] shows the following. We 
divided the [5th] density distribution in equipopulated bins of 
roughly 300 galaxies (x axis), and for each of these bins we plot 
the fraction of red galaxies in the first and last quartile of the 
[5g] density distribution (thick symbols are for the highest den- 



sity bin, light symbols for the lowest density bin). No trend with 
density on large scale is left, if we fix the small scale density. 

The bottom-right panel shows the opposite analysis. For any 
equipopulated bin of [5g] density {x axis), we plot the fraction 
of red galaxies in the first and last quartile of the [5th] density 
distribution in that bin. In this case, we see that a quite strong 
trend with density on small scale is left, if we fix the large scale 
density. This is more evident where [5th] reaches smaller scales, 
that is for higher densities. 

From these plots, we conclude that galaxy U - B colour is 
affected by small scale (< 1 /?"' Mpc) environment, while it ap- 
parently depends on the local density on larger scale only as a 
miiTor of the dependence on small scale. This at least it is true 
up to z = 0.5, and for a sample of galaxies with Mb < -18.5-z. 
Unfortunately it is not possible with our data to extend this study 
above this redshift, as we do not reach scales small enough (<1 
/j ' Mpc) given the increase of the mean interparticle separation. 



8. Comparison with other works 

Environmental effects on galaxy properties have already been 
studied within other data sets up to z ~ 1 and above. Considering 
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in particular spectroscopic surveys, the colour-density rela- 
tion has been analyzed within the VIMOS-VLT Deep Survey 
(VVDS) and the DEEP2 Galaxy Redshift Survey. Our results are 
qualitatively in agreement with their findings, showing that the 
colour-density relation in a luminosity selected sample weakens 
for increasing redshift. 

More in detail, using VV DS data (a pu rely flux limited 
galaxy sample with Iab < 24.0) Cucciati et al.! (2006) show that 
galaxy colour depends on the density contrast, with increasing 
strength for lower redshifts and for brighter galaxies. They show 
also the absence of the colour-density relation at z ~ 1, and a 
possible reversal of the relation itself for z > 1.2, for luminosi- 
ties around Mb(Ho=70) ^ -21.27. In the luminosity and redshift 
ranges that we have in common, we find that our results (F ig. [Hi 
are consistent with those presented in lCucciati et alJ(l2006i) . even 
if we are using here a slightly different density estimator. Our re- 
sults of Fig.|4]are in agree ment also with the on es obtained with 
DEEP2 data, as shown by .Cooper et alJ (l2007h when they use a 
purely luminosity selected sample (although their definition of 
red galaxy does depend on absolute magnitude). 

We remark that Cooper etal. (2007) obtain quite different 
results if they use a second differently selected sample. They 
show that the colour-density relation is still present at z ~ 1 
and disappears only at z ~ 1.3 if they use a luminosity-selected 
sample with a colour-dependent luminosity limit, the limit be- 
ing brighter for redder galaxies. In Section 16.11 we discussed 
the risk of analyzing the colour-density relation in a luminos- 
ity sele cted sample, a nd the colour-dependent luminosity limit 
used bv lCooper et al.l ([2007) may further decrease the control on 
biasing effects. The selection that they use includes only bluer 
galaxies at faint luminosities, while redder galaxies are included 
if they are bright. This means that the underlying spread in mass 
versus luminosity becomes luminosity (or colour) dependent, 
producing a non trivial eff'ect caused by the mass-density rela- 
tion. 

Our results on the EW [Oil] -density relation (Fig. |5] bot- 
tom panels) are also in agreeme nt with those obtained with the 
DEEP2 data ( ICooper et alj|2006l) . although we find a less signif- 
icant ^^[Olll-density relation within 0.85 < z < 1. 

Our findings are qu a litative ly in agreement also with those 
found bv lScovilleetal.1 (l2007al) . based on the entire COSMOS 
data set, with the density field computed using photometric red- 
shifts. In particular, they find that more massive objects reside 
preferentially in higher densities, where also the fraction of 
'Early Type' galaxies is higher, and this is true at all redshifts 
explored (0.2 < z < 1.1). Moreover, they find that the median 
value of TsF, defined as the inverse of the specific star forma- 
tion rate (sSFR), is higher in more dense environment, and this 
dependence weakens for higher redshifts. This is in agreement 
with our results on £^[011], that is a proxy for the sSFR (thus 

~ I/tsf). 

Finally, it is interesting to compare our results about the 
colour-density relation at fixed mass with similar works. In 
the local universe, Baldrv et al. (2006) find that the fraction of 
galaxies in the red sequence is higher in more dense environ- 
ments at any given mass in the range 9 < log(M/M0) < 11. 
[Kauffmann et al. (2004) show similar results considering the de- 
pendence of D4000A break and specific star form ation rate on 
local d ensity at fixed stellar mass. More in details, Baldrv et al.l 
(|2006|) find that the difference between the fraction of red se- 
quence galaxies in the lowest and highest densities decreases 
continuously for higher masses, as we also show in Fig.[TT] We 
note that their definition of "red sequence galaxies" implies a 
colour threshold that becomes redder for brighter luminosities. 



i.e. similar to the colour threshold we adopt, that becomes redder 
for higher masses. 

Considering higher redshift studies, we compared our re- 
sults on the colou r-mass-density relation to the ones presented in 
IScodeggio eFaP (I2009.) . based on VVDS data. While they find 
a mass-density relation for mass limited samples, they do not 
find any residual colour-density relation once fixing the stellar 
mass. In the redshift range that we have in common whit them 
(0.2 < z < 0.7) we can explore only stellar masses higher than 
the ones they study. Nevertheless, as we find a colour-density 
relation at fixed mass for 10.2 < log(M/M0) < 10.7 within 
0.1 < z < 0.5, while they do not find it for slightly lower stel- 
lar masses in the same redshift range, we studied in detail the 
differences between the two analyses. 

First, we observe that their red sequence in the colour-mass 
plane is much less populated then ours, with respect to the blue 
cloud. We verified this inspecting directly VV DS data, while 
for re ference we refer the reader to Fig. 1 in iFranzetti et al.l 
(l2007h . where a colour-magnitude diagram is shown, mirror- 
ing the colour-mass plane. De-populating randomly our red se- 
quence in order to lower our global fraction of red galaxies down 
to the value found in the VVDS sample, we find that the fraction 
of the remaining red galaxies is less depende nt on environment 
than before. Secondlv. lScodeggio et alJ(l2009l) use a density con- 
trast smoothed on a larger scale with respect to the one we use. 
Recomputing our colour-density relation in mass bins using the 
same scale they use, we obtain a weaker col our-density rela- 
tion. W e conclude that the differences between IScodeggio et al.l 
(l2009h and our work can be fairly explained, and that our results 
are robust against these differences. 



9. Discussion 

9.1. The evolution of the colour-density relation in a 
luminosity-selected sample 

It is nowadays well assessed that in the local universe clear con- 
nections exist between the local environment and several galaxy 
properties. On average, while redder and brighter galaxies live 
preferentially in regions where the local density is higher, the 



opposite is true for bluer and f a inter galaxies (e.;^. . i Balog h et aL 
2004; Ka uffmann etsll l2004t iBlanton et al.1 l2005l JHogg et aL 
2004). This picture is enriched by measurements of the two- 
point galaxy coiTelation function for different luminosity and 
colour sub-populations. More luminous galaxies tend to be more 
clustered that fainter ones, and red galaxies exhibit a stronger 
and steeper real-space correlation function than blue galaxie s 
(Norberg et al. 2002; Zehavi et al. 2002;'Madg wick et al.ll200l . 
Moreover, studying the luminosity function for early/red galax- 
ies in different environments, a brighter M* has been found for 
early-type galaxies residing in higher density regions, wh en den - 
sity is computed on several Mpc scales (Crotonetal. 2005h : 
note, however, that it is not yet clear whether the early-type 
LF computed in clusters has a brighter M* when compared 
with the LF in the fiel d (see for example the discussion in 
iBoselli & Gavazzill2006h . 

This composite picture at low redshift is globally still valid 
at higher redshift (z < 1). On the one hand, it has been 
found that the colour-density relation persists at higher red- 
shifts dCucciatietalJ 120061; ICooperet al1l2007h . On the other 
hand, studies of galaxy clustering by luminosity, b y spectral 
type and by stellar mass at hi gh redshift (e . ^. , Coil et al.l 
2004a b; Polio et al. 2006; Men eux et all l2006ll2008l but see 
also iMeneux et aLi2009i) have shown that red, early type galax- 
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ies at z ~ 1 are more strongly clustered than their blue counter- 
parts (this is also true at all redshifts up to z ~ 1 .2). Furthermore, 
the clustering length increases for more massive galaxies and for 
bright er ones. Thi s is m irrored by the LF per environment stud- 
ied in lllbert et al] (I200Q) . where it is shown that at z ~ 1 the 
brighter galaxies are preferentially populating high-density re- 
gions, while the opposite is seen for fainter galaxies. 

Although it is commonly observed that the pictures at z ~ 0.1 
and at z ~ 1 are linked by the progressive weakening of the en- 
vironmental e ffects for increasing redshift (ICucciati et al.ll2006t 
ICooper et al.l l2007). it is not still totally clear when this environ- 
mental dependence was established. It has even been observed 
the presence at z > 1 of colour-density and SFR-density rela- 
tions that are the opposite wit h respect to the oii e s in the local 
unive rse (ICucciati et al.1 l2006t lElbaz et al.1 120071: ICooper et al"] 
l2008h . Moreover, more light has still to be shed on the origin it- 
self of these environmental effects, to understand if galaxy prop- 
erties are related only to the environment in which galaxies have 
been born ('nature' hypothesis), or if the environment suiTound- 
ing galaxies at each epoch of their evolution continuously plays 
a role in shaping their properties ('nurture' hypothesis). 

Unfortunately, the evolution of the environmental effects on 
galaxy properties up to z ~ 1 and above has always been stud- 
ied in luminosity selected samples. We have shown in the pre- 
vious Sections that this selection may give rise to a misleading 
interpretation, given the underlying mass-density relation. For 
example, the faster decreasing with redshift of the red fraction in 
high d ensiti es than in low densities, found in both lCucciati et al] 
(l2006l) and ICooper et all (|2007), at a first level can be inter- 
preted as the proof that 'nurture' is at work, differently shap- 
ing galaxy properties in different environments. But when we 
take into account the mass-density relation, these findings be- 
come easily interpretable with a biased galaxy formation: i) 
mor e massive galaxies f ormed first in the highest density peaks 
(e.g. iMarinoni et al.ll2005h . and being more massive they con- 
sumed their gas reservoir more quickly, causing the fast increase 
of the red fraction with cosmic time in high densities; ii) low 
mass galaxies formed later in lower density regions, and we 
know that they consume their gas more slowly than their higher 
mass counterparts, and this is observable in the slow (or even 
absent) changing with time of the red fraction in low densities. 

Our results are qualitatively in agreement with the 
above-mentioned picture, when analysis is performed within 
luminosity-selected samples. In Sect. |4] for a luminosity- 
selected volume-limited sample complete up to z = 1, we find 
that the colour-density relation weakens for increasing redshift 
(Fig. lU, this trend being supported by the D„4000-density and 
fiW [Oil] -density relations (Fig. [Sj, that are clearly stronger in 
the lowest redshift bin explored (0.48 < z < 0.7) than in the 
highest one (0.85 < z < 1.0) 

More in detailed, we also find that the fraction of red galax- 
ies (but also of both 'passive' and 'active', as defined for Fig. 
|5]l in the highest density bin has a stronger variation with red- 
shift than the coiTesponding fraction in the lowest density bin. 
The above mentioned discussion about the galaxy biased for- 
mation and the mass-density relation is well suited to explain 
our findings. Nevertheless, it may be the case that the faster gas 
consumption in high densities is not only caused by the higher 
masses, but it may be accelerated by physical processes taking 
place in high density regions. It is not possible, within our re- 
sults, to distinguish which are the physical processes really at 
work in the overdense regions among the ones currently pro- 
posed (e.g. ram pressu re stripping of gas. JGunn & Gottlll972i 
galaxy-galaxy merging, iToomre & Toomral 19721 strangulation. 



iLarson et al.ll980l and harassment. lMoore et al.ll996l and so on). 
This prevents us from a more detailed analysis about the possible 
presence of nurture effects. 

Finally, we note that the redshift evolution that we find for 
the colour-density relation in a luminosity-selected sample is 
consiste nt with other paralle l analyses performed within the 10k- 
sample. llovino et aP (12010) show that the fraction of blue galax- 
ies Fhiue in galaxy groups is lower than Ft,iue among isolated 
galaxies at any redshift in the range 0.2 < z < 1.0, with this 
difference decreasing for higher redshift. Moreover, the fraction 
of morphologically early type is always higher in higher densi- 
ties, with environm ent parametrized u sing both the density field 
(iTasca et alJl2009l) and galaxy groups (iKovac etal .11201 Obh . 



9.2. The role of stellar mass in the colour-density relation 

Several studies indicate that, among other propert ies such as lu- 
minosity and morphological parameters, colour (Blanton et alj 
2005) and specific star formation rate (Kauffmanne t al. 20041) 
are the properties most tightly related to environment. In the past 
years it has been also shown that the gal axy stellar mass drives 
the star formation history of galaxies (e. e. iGavazzi & Scodeggid 
Il996t IScodeggio et al. 2002), that in turn is mainly responsible 
for the galaxy colour Moreover, the stellar mass is related o n 
the one side to the galaxy halo mass (iMandelbaum et aDl2006l) . 
and on the other side to the local environment as determined 
by nearby galax ies (Kauffmannet al. 2004; Bundy et al. 200^ 
IScod eggio et"an i2009l: iBolzonella et al.ll2009h. 

In particular, we refer the reader to^B olzonella et al] (l2009l) 
for a detailed study of the environmental effects on galaxy stel- 
lar mass within the zCOSMOS lOk-sample. They show that the 
Galaxy Stellar mass Functions (GSMF) is dominated by more 
massive galaxies in high densities, and that this difference be- 
tween the GSMF in lower and high density environments de- 
creases for higher redshift. We remark that these findings are 
miiToredby those related to the Luminosity Function (LF), based 
on the same data Zucca et al. (2009) . 

It follows that a non-biased analysis of environmental effects 
on galaxy colour (or on other properties) can be addressed only 
within complete mass-selected subsamples that span a narrow 
mass range. 

In our work, we disentangled the three-fold dependence 
among galaxy colour, mass and environment. Globally, we do 
not find any colour-density relation up to z ~ 1 when mass is 
fixed (Fig.[TO]l, with the exception of log(M/M0) < 10.7, where 
we see a clear colour-density relation (see also Fig.fTTI). As dis- 
cussed in Section[8] this is in broad agreement with simila r stud- 
ies in the local universe (Kauffmannet al. 2004; Baldr vetal] 
120061) . 

Our findings about the colour-density relation at fixed mass 
are coherent with the general picture depicted by other paral- 
lel analyses based on the zCOSMOS lOk-sample. iTasca et aP 
(2009) and Kovac et al. (2010b) find the same result studying the 
fraction of morphologically early type galaxies as a function of 
the density contrast and galaxy group membership, respectively. 
The same hold for the study of the fractio n of blue galaxies in 
groups and in isolation dlovino et al.ll2010l) . 

From these results, it seems that in the redshift range 0. 1 < 
z < 0.5 there exists a sort of threshold mass (log(M/MQ) ~ 10.7), 
below which the mix of the two broad galaxy populations (red 
and blue, early and late) does still depend on local density once 
stellar mass is fixed. Thus for low masses environment affects 
directly galaxy properties like at least colour and morphology. 
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Accordingly to these results, we asses that the colour de- 
pends primarily on mass, but for the low-mass regime the local 
environment modulates this dependence. 

First, we show that a colour-mass relation holds irrespec- 
tively of environment. Our definition of 'very red' galaxies im- 
plies a selection of redder galaxies for increasing mass. The fact 
that the fraction of 'very red' galaxies in high densities is al- 
most constant as a function of mass (thick circles in the mid- 
dle panel of Fig. fTTT i means that in this environment galaxy 
colour becomes globally redder for higher masses. This effect 
is even stronger in low density regions, where the 'very red' 
ga laxy faction increases f or increasing mass. The same is shown 
in iBolzonella et al.ll2009l studying the GSMF per galaxy types 
in different environments. In parall el, this behaviour is also ob- 
served in the Luminosity Function dZucca et al.ll2009n . 

Second, in this scenario where a colour-mass relation is em- 
bedded in both low and high density regions, the environment 
seems to have a role in shaping different colour distributions in 
the redshift range 0.1 < z < 0.5 once mass is fixed. This means 
that environment does play a role independently of stellar mass. 
For example, in the central panel of Fig. [TT] we see that the frac- 
tion of red galaxies is clearly a function of the density contrast at 
any given mass below log(M/M0) ~ 10.7 (the possible presence 
of a similar trend for log(M/Mo) ~ 11 . 1 is discussed in Sect. 
I6.2I 1. Continuin g the comparison with th e parallel analysis on the 
GSMF, also in IBolzonella et alj (12009 ) the role of environment 
on galaxy colour per mass bin emerges clearly. The fractional 
contribution of the separated blue and red galaxies GSMFs to the 
global GSMF varies with environment, especially in the range of 
intermediate masses (9.5 < log(M/Mo) < 10.5). 

A possible general picture to describe these findings is the 
following. Considering that the reddest and most massive galax- 
ies have been formed earlier (at z > 2, having akeady an age 
of > 1 Gyr at z ~ 1, see e.g. Vergani et al. 2008] and references 
therein), and that moreover their Star Formation Histories ( SFH) 
are typically faster (< 3 Gyr, see e.g. iGavazzi et al.ll2002h than 
those of lower mass galaxies, their formation epoch and their 
evolution time-scale took place on the mean before an appre- 
ciable growth of structures (at z ~ 1 the number of structures 
with log(M/M0) > 5 X 10'"* was ~ 100 times lowe r than today, 
in the concordance cosmology, see lBorganill2006l) . We do not 
see evident environmental effects on these galaxies because en- 
vironment could not affect their evolution. On the contrary, lower 
mass galaxies (log(M/Mo) < 10.7) not only have been formed 
more recently, but their SFH are also slower Therefore, physical 
processes typical of high density regions could have modified 
their SFH itself. For a broad picture about the role of mass and 
environment in driving galaxy evolution, obtained w i th SP SS 
and zCOSMOS data, we refer the reader to lPeng et all (l2010h . 

This picture is in agreement with the scenario that has been 
already proposed in literature where bot h 'nature' and 'nur- 
ture' affect galaxy p r operties (see e .g. Kauff mann et al.l 12004 
iDe Lucia et al.1 l2006l ICucciati et al.l |2006.) . to which we can 
add t he dependence on m ass. First, galaxy formation is 'bi- 
ased' (iMarinoni et al.l l2005'). meaning that more massive galax- 
ies formed first in the highest density peaks, later on followed 
by lower mass galaxies in less dense environment. Second, for 
relatively low mass galaxies the evolution is affected by com- 
plex physical processes that depend on the local environment. 
According to this picture, the properties of a low-mass galaxy 
are thus affected by the environment in which the galaxy resides 
at any given epoch of its evolution, and not only at the time of 
its formation, as a simple imprinting. 



Unfortunately, the relatively large uncertainties in our mea- 
surements prevent us from assessing in a precise and quantitative 
way the relative roles that 'nature' and 'nurture' did play. We de- 
fer this study to a future work. 



10. Summary and conclusions 

In this work we use the first ~ 10000 spectra of the purely 
flux limited zCOSMOS bright sample (Iab < 22.5) to inves- 
tigate the redshift evolution of environmental effects on galaxy 
spectro-photometric properties, mainly U-B colour, in the range 
0.1 < z < 1.0. We refer the reader to Kovac et al. (2010a) for a 
full description of the environment parametrization via the local 
density field. More in details, we want to disentangle the mul- 
tiple dependence among local density, galaxy colour and stellar 
mass, in order to verify whether the environment acts directly on 
both stellar mass and colour, that are related to each other. Our 
main results are the following. 

- Using a luminosity-selected volume-limited sample (Mb < 
-20.5 - z) we find, confirming previous results, that the frac- 
tion of red (U-B > 1) galaxies depends on environment at 
least up to z ~ 1, with red galaxies residing preferentially in 
high density environments. This trend becomes weaker for 
higher redshifts, and it is mirrored by the one we observe for 
the fraction of galaxies with D„4000 > 1 .4. We also find that 
the fraction of galaxies with log(EW[OII]) > 1.15 is higher 
for lower densities, and this is true up to z ~ 1 . Also in this 
case we observe a weakening of the environmental effects 
for increasing redshift. 

- Giv en the fact that ste llar mass depends on local density (see 
e.g. I Bolzonella et alJ 1200 9 and our Fig. |7]i, and given that 
colour depends on mass, it is difficult to interpret the mean- 
ing of the colour-density relation in a luminosity-selected 
sample. In fact, the wide spread in mass-to-light ratios em- 
bedded in this selection leave us with a broad range in stellar 
masses, that biases any possible direct colour-density rela- 
tion via the mass segregation as a function of environment. 

- We disentangle the colour-mass-density relation using 
galaxy subsamples enclosed in narrow mass bins, to avoid 
any dependence of mass on local density. We study the 
colour-density relation in mass bins of A log(M/M0) = 0.2, 
in three different redshift bins within the range 0.1-1.0, pay- 
ing attention to attain to the mass limits imposed by the flux 
limit of our survey. We find that , once the mass is fixed, the 
colour-density relation is globally flat, but with some excep- 
tions. We observe that within 0.1 < z < 0.5 the fraction of 
'very red' galaxies depends on environment even when mass 
is fixed, at least for log(M/M0) < 10.7 (see Section |6^ for 
the definition of 'very red' galaxies). This means that en- 
vironment affects directly not only the stellar mass, but also 
other galaxy properties, at least for these given mass and red- 
shift ranges. 

According to these results, we suggest a scenario in which 
the colour depends primarily on mass, but for a relatively low- 
mass regime (10.2 < log(M/M0) < 10.7) the local environment 
modulates this dependence. These galaxies have been formed 
more recently, in an epoch where evolved structures were al- 
ready in place. Moreover, the physical processes typical of high 
density regions could operate during longer periods of time in 
shaping the properties of these galaxies, thanks to their longer 
(on average) star formation histories. 
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It is not possible with this analysis to quantitatively assess 
the relative role of biased initial conditions (more massive galax- 
ies form first and in more dense environments) and of physical 
processes acting during galaxies lifetimes, in order to precisely 
explain the origin (or the absence) of both colour- and mass- 
density relations. We refer the reader to lPenget alj (1201 Oh for a 
more comprehensive analysis on how mass and environment af- 
fect galaxy evolution. On the one hand, the forthcoming zCOS- 
MOS 20k bright survey will allow us a more detailed analysis of 
these mass and redshift regimes, on the other hand this work has 
to be complemented by the study of the redshift evolution of the 
less massive galaxies population, that is not possible within the 
zCOSMOS bright sample, and that is of crucial importance to 
better understand how environment affects the different galaxy 
populations. 
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